Cascading Gravity and Degravitation 

Claudia de Rham 1,2 , Stefan Hofmann 1,3 , Justin Khoury 1 , Andrew J. Tolley 1 

1 Perimeter Institute for Theoretical Physics, 
31 Caroline St. N., Waterloo, ON, N2L 2Y5, Canada 
2 Dept. of Physics & Astronomy, McMaster University, Hamilton ON, Canada 
3 NORDITA, Roslagstullsbacken 23, 106 91 Stockholm, Sweden 

Abstract 



We construct a cascading brane model of gravity in which the behavior of the gravitational 
force law interpolates from (n + 4)-dimensional to (n + 3)-dimcnsional all the way down to 
4-dimensional from longer to shorter length scales. We show that at the linearized level, 
this model exhibits the features necessary for degravitation of the cosmological constant. 
The model is shown to be ghost free with the addition of suitable brane kinetic operators, 
and we demonstrate this using a number of independent procedures. Consequently this is 
a consistent IR modification of gravity, providing a promising framework for a dynamical, 
degravitating solution of the cosmological constant problem. 
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I. INTRODUCTION 

The remarkable theoretical challenge posed by the cosmological constant problem [l| has spurred 
many attempts at directly modifying Einstein's gravity at large distances. A compelling example 
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of one such infrared (IR) modification is the DGP brane- world model ln this scenario, our 

visible world is confined to a brane in an infinite 5D bulk. The inclusion of a graviton kinetic term 
on the brane recovers the usual gravitational force law scaling, 1/r 2 , at short distances, but at 
large distances it asymptotes to the 5D scaling, 1/r 3 . 

The resulting weakening of gravity on large scales can have profound implications for the cos- 
mological constant problem Q|. In particular, these models have in common the feature that the 
4D graviton acquires a mass, or is a resonance — an infinite superposition of massive states. In a 
recent proposal [jj], two of us put forward the idea, based on earlier work in that if gravity is 
sufficiently weakened in the infrared, then vacuum energy could effectively decouple from gravity 
or degravitate over time. This would provide a causal, dynamical solution to the cosmological con- 
stant problem. However, the degravitation phenomenon is not realized in the standard DGP model 
because the weakening of gravity is not sufficiently steep. In [5j it was shown that degravitation 
requires that the AD graviton propagator grows slower than 1/p in the IR [5], whereas in standard 
DGP it grows exactly as 1/p. An explicit model with the desired behavior for the propagator was 
not given in [5|, and the goal of the present work is to construct specific non-linear theories that 
can exhibit the degravitation phenomenon. 

A promising way to realize degravitation is to consider higher-codimension DGP brane-world 
models, i.e. higher-dimensional bulks. If the bulk has six or more space-time dimensions, then 
the gravitational force law falls off faster in the IR than in the standard DGP model, and as we 
shall see in what follows this fall off is sufficiently fast to exhibit degravitation at the linearized 
level. However, higher-codimension DGP scenarios have proven notoriously difficult to realize 
consistently for two reasons: Firstly, the 4D propagator diverges when evaluated on the brane 
and must be regularised flfl Q. Secondly, the simplest nonactions are plagned by ghost-like 
instabilities [101, LLH • 

In this paper we argue that both pathologies — divergent propagator and ghost instabilities 
- are resolved in the recently proposed Cascading DGP model 12|]. The idea is that our visible 
3-brane lies within a succession of higher-dimensional branes, each with their own induced grav- 
ity terms, embedded in one another in a flat D-dimensional bulk. For instance, in the simplest 
codimension-two case, our 3-brane is embedded in a 4-brane within a flat 6D bulk. The gravita- 
tional force law therefore "cascades" from \D (1/r 2 ) to 5D (1/r 3 ) to GD (1/r 4 ) etc., as we probe 
larger distances on the 3-brane (a similar cascading behavior of the force law was also obtained 
recently in a different codimension-two framework 
branes was discussed in 1J| with somewhat different motivations). 



131 ]. and closely related work on intersecting 



In Sec. [H] we review the idea of degravitation to motivate the study of higher-codimension 
brane set-ups, and explain how these satisfy the properties expected of the linear theory using 
the spectral representation of the propagator, namely that the 4D propagator in momentum space 
grows more slowly than 1/p as p — > 0. Cascading DGP gravity, therefore, provides a candidate 
nonlinear realization of the degravitation idea and opens the door to a dynamical resolution to the 
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cosmological constant problem. 

We begin with a scalar analogue in Sec. IIIH consisting of a bulk scalar field with induced kinetic 
terms on each of the branes. We find that the usual divergence of the 4D propagator is cured by 
the presence of the higher-dimensional branes — the extra kinetic terms act as regulators for the 
3-brane propagator. This framework is then generalized to gravity in Sec. IIV| where we find that 
the cascading set-up similarly regulates the usual divergences, however it does not automatically 
avoid the ghost instability characteristic of higher-codimension DGP models [lol . [ill ]. As in the 
pure codimension-n DGP model, adding an R4 term on the 3-brane results in a ghost instability. 
This can readily be seen at the level of the one-particle exchange amplitude, which naturally splits 
into a transverse-traceless (or massive graviton) contribution and a scalar component. In the IR, 
these combine to give the usual tensor structure of a massless graviton in D dimensions, as it 
should. In the UV, however, the scalar contribution is negative, signaling a ghost, i.e. the scalar 
part of the Green's function does not allow a well-defined spectral representation. 

In 12] it was argued in the co dimension- two case that the ghost is cured by including a suffi- 
ciently large tension on the (flat) 3-brane. Here we follow a different approach by noting that an 
i?4 term is not the most general induced gravity action on a 3-brane. In Sec. [V] we present two 
alternatives, obtained through regularization of the codimension-n brane. An important concern 
when smearing out branes is maintaining gauge invariance, an issue we take great care in by util- 
ising the Stiickelberg formalism. The approach of Sec. |V]is not restricted to Cascading DGP, but 
applies more generally to arbitrary codimension-n DGP models. 

• The first method (Sec. IV Bp consists of blowing up the 3-brane into a codimension-one sphere 
of radius A. The key observation is that, once the defect has finite thickness, its worldvolume 
has more than 3+1 dimensions. Consequently, it no longer makes sense to include an R4 
term covariantly. Rather, the natural choice is to add an R3+ n term, the Einstein-Hilbert 
term appropriate for the dimensionality of the regularized codimension-one brane. 

• The second method (Sec. IV Cj) promotes the 3-brane into a full (3 + n)-dimensional object of 
radius A. The natural choice in this case is to include an Rd term within its worldvolume, 
albeit with a different Newton's constant than that in the bulk. Thus the codimension-n is 



understood as a medium with non-zero gravitational permeability, as first introduced in 151 ] . 



Remarkably, these additional kinetic terms on the regularized brane have dramatic consequences 
for the scalar degrees of freedom: Both methods yield the same, ghost-free exchange amplitude, 
displaying the tensor structure characteristic of the D-dimensional massless graviton on all scales. 
The effect of these terms is very similar to those considered in [ll|] where it was shown that adding 
an Rq kinetic term on a 3-brane in 6D would, for sufficiently large value of this term, remove the 
ghost. Another way to understand these operators is that one is adding extrinsic curvature terms 
into the codimension-two brane action. While these would normally render the boundary value 
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problem ill-defined without the addition of higher derivative terms in the bulk, here we can make 
sense of these terms by suitably regularizing the brane. 

That these two methods yield the same answer can be understood most transparently in the 
language of effective field theory, as shown in Sec. IV Dl In the thin-brane limit, A — > 0, we can 
perform a Kaluza-Klein reduction within the worldvolume of our defect to obtain an effective action 
on the 3-brane. Note that dimensional reduction is performed only within the defect, while the 
bulk remains infinite in extent. The resulting effective action reduces to an R<± term, although with 
additional scalar-tensor couplings to the extra-dimensional moduli. These scalar-tensor couplings 
are essential in modifying the scalar part of the 4D exchange amplitude, thereby rendering it 
ghost-free. 

As mentioned earlier, these regularization schemes apply quite generally to any higher- 
codimension DGP model. In Sec. IVI} however, we focus on their implications for cascading DGP. 
In the codimension-two case, in particular, the one-particle exchange amplitude on the 3-brane 
interpolates between the usual 1/r at short distances and 1/r 3 at large distances, while the tensor 
structure remains that of a 6D massless graviton on all scales. 

The decoupling limit 1(| of the codimension-two cascading theory offers yet another perspective 
on how the ghost is excised from the spectrum, as shown in Sec. IVIII By taking M§,Mq — > oo 

1 /7 

keeping the strong coupling scale A s = (Mg 6 /M|) fixed, we can integrate out the 6D bulk 
physics to obtain a local, five-dimensional effective theory for the 5D metric and the scalar (brane- 
bending) mode tt [12]. For the usual choice of R4 term on the 3-brane, the tt scalar ends up having 
a negative kinetic term on the 3-brane — hence it is a ghost. On the other hand, the 3-brane 
effective action mentioned above contributes a positive induced kinetic term for tt, rendering this 
mode healthy. In fact, the decoupling analysis shows that our effective action is part of a broader 



class of 3-brane lagrangians invo 
non-zero tension considered in 



ving tt that yield a ghost-free amplitude. The inclusion of a 
12| is another example. 

At face value the 6D tensor structure of the 4Z) graviton propagator is inconsistent with solar 
system tests of General Relativity. In other words, the codimension-two Cascading DGP framework 
yields —1/4 as the TT' coefficient in the exchange amplitude, whereas General Relativity yields 
— 1/2. This mismatch persists for arbitrarily small graviton mass, a paradox known as the van 
Dam-Veltman-Zakharov (vDVZ) discontinuity [l3]. As in the standard DGP model, we expect 
this effect to be an artifact of perturbation theory [lij. In the vicinity of astrophysical sources, 
within the so-called Vainshtein radius, the one-particle exchange approximation breaks down and 
the scalar modes become strongly coupled. The end result is that the scalar modes decouple, and 
Einstein gravity is a good approximation. 

As argued in Sec. IVIIIl the natural outcome in Cascading DGP is a succession of Vainstein 
radii for widely separated cross-over scales. In the codimension-two case, for instance, the scalar 
responsible for the 6D — ► 5D transition is tt. We therefore expect this scalar to decouple below the 
Vainshtein radius as in the standard DGP model in 6D. Below this scale, the theory effectively 
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reduces to 5D gravity coupled to a 3-brane with induced gravity term, i.e., the standard DGP 
model. We can anticipate a second Vainshtein effect for the longitudinal mode of the graviton at 
sufficiently short distances. 

To summarize, in this paper we show that the cascading DGP framework recently proposed 
in 12| is a compelling candidate for a consistent nonlinear realization of degravitation of the 
cosmological constant. Motivated by independent regularization schemes, we derive a new effective 
action on the 3-brane, which eradicates the ghost instability characteristic of higher codimension 
DGP models. In the codimension- two case, the resulting gravitational force law on the brane 
interpolates from 6D at the largest scales to 4D at the shortest via an intermediate 5D regime. 
This model naturally generalizes to higher dimensions where there is a similar cascade from (4+n)- 
D to (3 + n)-D down to 4D. 

There are many important issues which are left for future work. First of all, it remains to 
show that degravitation takes place in the full non-linear theory, by studying the cosmological 
evolution in the presence of large vacuum energy on the brane. The linearized analysis presented 
here certainly displays all the desired features, but it remains to be seen whether this is also the 
case at the non-linear level. Moreover, one should check in detail that the extra scalar modes 
decouple near large sources due to non- linear (or Vainshtein) effects, a necessary condition for the 
theory to be phenomenologically viable. 

II. REVIEW OF THE DEGRAVITATION PHENOMENON 

Almost all attempts in addressing the cosmological constant problem have focused on making 
the vacuum energy density small. Degravitation is a framework that instead allows for a large 
cosmological constant, but suppresses its backreaction by making gravity exponentially weaker at 
long distances. In other words, the vacuum energy density is not small in this picture, it is simply 
degravitated. A phenomenological modification to Einstein gravity that encapsulates this behavior 

Q 



is 



G^ x {L 2 U)G 11v = %ttT 11v , (1) 



where Newton's constant has been promoted to a derivative operator. This is such that gravity 
behaves as a high-pass filter with characteristic scale L: Gn — > G-S for L 2 (— □) — ► oo and — > for 
L 2 (— □) — ► 0. However, as it stands this equation cannot be consistent. This is immediately seen, 
for instance, by noting that the Bianchi identity — a direct consequence of general covariance in 
four dimensions — is violated. 

Indeed, a key observation made in [j], is that any degravitating theory must reduce at the 
linearized level to a theory in which the graviton has a mass or is a resonance, i.e. a continuum of 
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massive states. That is, the linearized equations of motion must be of the form 

(^U + (V - V^h) = —T^v , (2) 

where (£h) = —\Dhn U /2 + . . . is the linearized Einstein tensor. The above is a straightforward 



generalization of Fierz-Pauli massive gravity [201 ] . the only ghost-free theory of a (free) massive spin- 
2 particle, with the mass term promoted to a function of the derivative operator. An immediate 
corollary is that any degravitating theory must describe extra degrees of freedom, corresponding 
to the extra polarizations of the massive graviton. Indeed, the metric fluctuation hau propagates 
5 degrees of freedom: 2 helicity-2 modes, 2 helicity-1 modes, and 1 helicity-0 mode. 

Deconstructing the full metric fluctuation into its various helicity constituents offers a beautiful 
proof of concept for degravitation: the various helicities are honest gauge fields, i.e. they transform 
in a nontrivial way under small gauge shifts. However, the full metric fluctuation is gauge invariant. 
In other words, the helicity-1 and helicity-0 components play the role of Stiickelberg fields that 
render the full metric fluctuation gauge invariant. Thus, hau becomes an observable (at least 
in principle), and is an equation of motion for an observable quantity. Hence, a sensible 
solution of ([2]) has to be bounded. This is not the case for a massless fluctuation which, for a 
homogeneous source, develops an instability towards de Sitter. In contrast, massive fluctuations 
respect Minkowski, i.e. the mass term in ([2]) guarantees that flat space is a consistent background, 
even for a cosmological term. At this point we reach the conclusion that the mass term in ([2]) gives 
rise to a dispersion for the gravitons of generalized Fierz-Pauli that effectively decouples them from 
any infrared source. In other words, a high-pass filter must be at work and is activated under the 

aquation 

t Q, Q- 

The next question is what are the allowed forms for m 2 (D)? Since vacuum energy is the 
longest-wavelength source, we are interested in the far IR behavior of this mass operator. A useful 
parametrization in this regime is a power-law: 

m 2 {U) = L- 2 ^ l ~^U a . (3) 

Clearly a < 1 in order for this to represent an IR modification of gravity. Requiring that the 
corresponding propagator has a well-defined spectral representation, so that the theory is free of 
ghosts, puts a lower bound on this parameter: a > 0. Thus all degravitation theories can be 
classified at the linearized level by a single parameter a within the allowed range: 

< a < 1 . (4) 

For instance, massive gravity corresponds to a = 0, whereas the DGP model has a = 1/2, as can be 
seen from f)18|) for instance. Moreover, the logarithmic behavior as p — > of our codimension-two 
cascading propagator, as shown in (f25j) . essentially corresponds to a ps 0. 



spell of gauge invariance. Indeed, the filter equation (TjQ) was shown to arise as an effective equation 
for the helicity-2 (or Einsteinian) modes, once the other helicities have been integrated out 
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A. Spectral Representations 

We now present an argument that all cases in which the gravitational potential at large distances 
behaves as higher dimensional, correspond to a = or a = 1/2. Associated with a given spectral 
representation of the propagator, 

we have the 'Newtonian potential' 

V(r)= r p ( S ) eXp( -^V (6) 
Jo r 

In D space-time dimensions, the Newtonian potential scales as V(r) ~ l/r D ~ 3 , and so in a 
codimension-n DGP model, we expect that at large distances we feel the full (4 + n)-dimensions 
so that V(r) ~ l/r 1+n . Since the potential at large distances is determined by the s — ► behavior 
of the spectral density, it is clear that p(s) must scale as p(s) ~ W 2 " 1 as s — > 0. To make this 
concrete, define p(s) = s n ^ 2 ~ 1 f(s), where the only requirement is that the function f(s) is finite 
and nonzero at s = 0. The Newtonian potential is then 



Jo r 



which can be simplified by redefining the integration variable s = x jr so that 

2 roc 
r~ - jo 

On taking the limit r-> oo we have 



V(r) = I x n - l f{x 2 /r 2 ) exp (-x)dx . (8) 



as anticipated. 

We are now in a position to show that all higher codimension cases corresponds to a = 0. In 
the limit p — > 0, the propagator is given by 

poo s n/2-l f( \ 

G(0)= / J -^ds. (10) 

Jo s 

This integral is convergent for small s provided that n > 2, and hence 1/Q(p^) ~ p 2a with a = 
for all codimensions n > 2. For the case n = 2 the integral is IR log divergent, G(Pu) ~ l°g(p), 



which is consistent with the above results. For the case n = 1 we have 

r-OO ,-1/2 

Jn=\" 



By redefining s = p 2 y, we find 

r°° ,.-1/2 f f„2 \ /-oo -1/2 

lim g^x^) = limp- 1 / y -f^ dy = p" 1 /^) / ^— = p~ V(0), (12) 

It follows that n = 1 always corresponds to a = 1/2. This is familiar from the usual DGP result. 
We stress that all of these statements depend only on the leading behavior of the spectral density 
as s — > and so are independent of other details of the set-up. 
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III. CASCADING SCALARS 

In this Section we review the recently proposed cascading mechanism [12] to regularize physics 
on a codimension-two brane, and generalize it to arbitrary codimensions. For simplicity we focus 
here on a toy model of a scalar field, leaving the generalization to gravity to Sees. IIVI and IVI1 
The idea is to embed the codimension-two object on a codimension-one brane, each with their 
own intrinsic kinetic terms. (In the gravitational case, these are intrinsic Einstein-Hilbert terms, 
as in the DGP model.) The resulting force law between test particles on the codimension-two 
brane is therefore 4D at short distances, then effectively 5D over some range of scales, and finally 
6D at large distances. We also show how this mechanism is completely generalizable to arbitrary 
codimensions. 

In what follows, we use the notation that capital Latin indices run over all the D dimensions, 
A,B = 0, 1, . . . , (D — 1), small Latin indices run over the n = (D — 4) extra dimensions, a,b = 
4, . . . , (D — 1), while greek indices span over our four dimensions [a, v = 0, . . . 3. When working in 
D = 6 dimensions, we will use the additional notation that greek indices a, (5 span over the five 
dimensions of a codimension-one brane. designates the (i-dimensional "Planck" scale. 



A. Codimension-two case 



In the codimension-two case, the simplest construction is to consider a (3 + l)-brane embedded 
on a (4 + l)-brane, which itself is a codimension-one object in a 6D bulk. Denoting the 6th 
dimension coordinate by z, and the 5th by y, the (4 + l)-brane is taken to be at z = 0, while the 
(3 + l)-brane on its worldvolume is again at y = z = 0. A Z2 symmetry is assumed across both 
branes, so the y and z-coordinate range from to 00. Unlike the pure codimension-two case, here 
we will find that the propagator is finite when evaluated on the codimension-two brane. That is, the 
embedding on a codimension-one object removes the usual UV divergences of higher-codimension 
brane theories, 

aaa. 

The full action is given by 

S = — ^ d 6 xd A $d A <S> - -5- / d 5 xd a $d a <S> - -± j d^xd^® , (13) 

where the scalar field $ is dimensionless. The induced propagator between sources on the brane is 
usually calculated using the boundary effective action approach — see, e.g., [if]]. Here we follow a 
different path, by treating the brane action as a localized coupling for the bulk free theory 0, Q]: 

S = -^Ljd 6 x d A t>d A t> - j d e x + ^Sgtog} , (14) 



with 

A/f 3 M 2 



£^>nng = d a $d a $ ; ££&tag = -^5(z)S(y) d^P* . (15) 
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This corresponds to a coupling Ai = —M^(qf +p 2 ) on the codimension-one brane [z = 0), where q\ 
is the momentum in the y direction, p 2 = rj^PaPv and p^ the 4-momentum in the codimension-two 
brane directions x M , along with a coupling A2 = —M 2 p 2 localized on the codimension-two brane 
(z = y = 0). 

The free propagator associated with the bulk action is simply 

1 f dq e iq ( z ~ z "> 

Goip^qDz,^) = —j / — — — — — 2 ' ( 16 ) 

Ml J 7T p z + qf + q z 

where we work in the mixed representation, i.e., in real space along z and in momentum space 
along the y and directions, and q is the momentum conjugate to z. When both legs are taken 
on the codimension-one brane (z = z' = 0), this propagator reduces to 

£o(iv, qi) = Goip^ qr, 0,0) = f — 2 I 2 = 1 = . (17) 

Notice that in this scenario, due to the Z2 symmetry imposed across both the brane, only half of 
the real plane for z is considered. This explains the normalization taken for q. 

i) Coupling on the codimension-one brane 

This free propagator is modified by the interaction term on the codimension-one brane. These 
corrections are symbolically represented in Fig.[TJ where in n = 6 dimensions the coupling constant, 
Ai = -Ml(p 2 +ql), is read off from -Coupling- 

Summing up these diagrams, we obtain the modified 

Green's function Gi on the brane: 



Gi{pn,qi) = Go(Pn,qi) + \iGo(Pfi,qi) 2 (l + \iGo(Pn,qi) + (\iGo(p^,qi)) 

goQ^gi) = _L 1 

1 + Ml (p2 + qf) g {p^ qi ) Ml p 2 + q f + m^^T^qJ ' 

which is recognized as the usual DGP propagator, with 5D to 6D crossover scale 



2 + 



(18) 



Mi 

To proceed further it is useful to transform back to real space in y, so that 

Gi{p^y,y ) = ^rjs — 7. 7= 9 ■ (20) 

On the codimension-two brane, we therefore have 



Gi(p„) = 0,0) = < 



k ^^^ tm-1 (v*S) for p>m6 ' 



(21) 
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Notice that this integral is only finite as long as M5 is non-zero. Physically, the codimension- 
one brane plays the role of a regulator without which the Green's function on the co dimension- two 
brane would be ill-defined. In particular, in the limit M5 —* one recovers the logarithmic diver- 
gence, Qiipn) ~ J Q A d<li — 4/2 2 ~ ]°e{A-/p), characteristic of codimension-two branes 0]. Thus 
including a kinetic term on the codimension-one brane is crucial in order to avoid any singularity, 
making the scalar field well-defined on the codimension-two brane. In the gravitational case, the 
presence of kinetic terms on the codimension-two brane will also be sufficient to avoid a singularity, 
as we will see in Sec. IIVI 

ii) Coupling on the codimension-two brane 

The same technique can be applied once more to take into account the coupling ^coupling local- 
ized on the codimension-two brane. The diagrams contributing to the new Green's function are as 
shown in Fig.[U with coupling A2 = —M 2 p 2 localized at (y, z) = (0, 0). Using the same summation, 
we therefore obtain the final Green's function on the codimension-two brane 

a r„ \ - = 1 1 (99) 

1 + M 2 p 2 01 (P„) M| p 2 + 5(p 2) ' W 



where for future reference we have defined 



g(p 2 



irm 5 y/™%-P 2 



4 tanh -i ( /™sEP 



tan" 1 ! 



for p < me 



for p > ijiq , 



(23) 



with ins the 4D to 5D cross-over scale, 

m 5 = M . (24) 
5 M 2 y ' 

Once again, in the absence of the original codimension-one brane, corresponding to setting M5 = 

or, equivalently, m% 3> p, the propagator on the brane would be ill-defined: 

Q2 iPy) 7^4 log p/m 6 . (25) 

Thus the codimension-one brane plays the role of a regulator, allowing for an arbitrary source to 
be localized on the codimension-two brane. 

We can also notice that we recover the DGP limit as m% — ► 0, i.e. when the field $ effectively 
only propagates in five dimensions: 

p»me 1 1 

M| p 2 + m^p 

The previous example therefore exhibits a "cascading" behavior, interpolating from a four- 
dimensional behavior at short distances, p 3> ms, to a five-dimensional one for intermediate scales, 
m5 3> p S> m6, and finally a six-dimensional behavior at large distances, p <C mj. If > ijiq 
then there is a direct transition from 4D to 6D at a scale p ~ y/m^rriQ. 
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A m A m A m 

FIG. 1: Coupling corrections to the two-point function. The Green's function on the codimension-m brane 
is related to that on the codimension- (tti — 1), Gm—i(p)- A m is the coupling associated with -Coupling- 
A m = -MtXlZ™ (p 2 + q\ + ■ ■ ■ + ql_ m ) ■ 

B. General Codimension- n 

As can already been inferred at this point, this method can easily be generalized to any higher 
codimension. In particular for a codimension-n brane, the full action will be given by 
M n+2 f M n+3 f M 2 f 

s = — |ti / d n+A xd M $>d M § — |^ / d n+3 xd a $d a <s> / d A xd^d^ , (27) 

and a coupling A m = —M^™Zm(P 2 + Qi + ' ' " + Qn-m) wm be present on each codimension-m brane 
(0 < m < n — 1). Here is the momentum along the (3 + l)-brane and q m the momentum along 
the m th extra dimension. We denote by Go the free (4 + n)-dimensional Green's function 

GoiP^Qir ■ ■ ,q n -i) = , (28) 

MZtlyJp 2 + 9? + • • • + Q 2 n-i 

evaluated on the codimension-one brane Go = Go{z n = z' n = 0), and by Q m the Green's function 
evaluated on the codimension-m brane that takes into account the couplings Ai,--- ,A m . This 
Green's function is thus related to that on the codimension- (m — 1), G m ~i, by 



—m Gm—l 

ym{Pn,qi, ■ ■ ■ ,Qn-m-l) = / Z r -p, 7 

J ^ 1 - A m t/ m _i(p M ;gi,- ■ ■ ,q n - r 







f dq n - m Grn-liPfj,] qi, ■ ■ ■ , gn-m) 

~ J * 1 + M%+%I™(p 2 +$£ + •■■ + q 2 n .J Gm-1 (p„; qi, • ■ ■ , qn-m) ' 
It follows that the field propagator projected on the codimension-m (3 + l)-brane is simply given 
by 

*w = . +j &&Ur (30) 

with Gn-i satisfying the previous recursive relation. This gives rise to an arbitrary dimension cas- 
cading setup where the field progressively switches from a four-dimensional to a (4 + n)-dimensional 
behaviour as larger distances are considered. 

IV. CASCADING GRAVITY 



We would like to generalize this scalar field toy model to the gravitational case, where brane 
induced kinetic terms are now replaced by Einstein-Hilbert terms. In this section, we will focus 
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on the co dimension- two case, although the generalization to any codimension is straightforward. 
In the scalar field toy-model, the presence of a kinetic term on the codimension-one brane effec- 
tively regularizes the induced codimension- two propagator. The situation for gravity, while similar, 
contains some subtleties. 

In the six-dimensional scenario, the straightforward generalization of (|13j) to gravity is 

S = -6 J # x J=^Rs + d 5 xV^R5 + d 4 x^9~ A R A . (31) 

In this case, the tensor modes behave precisely as the scalar field of Sec. 1111} as expected. In 
particular, its induced propagator in four dimensions is made finite by the embedding on the 
codimension-one brane. However, a key difference, as we will see, is that one of the scalar modes 
propagates a ghost 

We first present in what follows the aforementioned pathology associated with (|3ip . The detailed 
derivation of the gravitational exchange amplitude is presented in Appendix^] We summarize here 
the main result, mainly focusing on the behavior of the scalar modes, since they are the ones that 
have a ghost-like behavior at high energies. 

Working around flat space-time, the perturbations may be decomposed into a scalar part and 
a four-dimensional transverse and traceless (TT) part, 

hfw = hjj + irr/^ + gauge terms . (32) 

The tensor modes only couple to the conserved and traceless part of the source, i.e. to = 
T$ — ^T^r]^ + |^T( 4 ), and their propagator is precisely that of the scalar field toy model of 
section Sec. IIIII On the codimension-one brane, their equation of motion is governed by 



Mf 



□ 5 - m 6 hJJ = 5{y) ( + -fU 4 h™ ) . (33) 



2 

Their propagator is given in (|22|) . and we obtain 

g = 26 2W S^ = ^ _ n4+ 1 g( _ D4) IV (34) 

Notice that since these modes behave identically to the scalar field toy-model, they are completely 
regularized by the codimension-one brane. 

The scalar mode, ir, is also regularized by the codimension-one brane, 



Mf r 



□ B - m 6 v/=Oi] tt = 1 5{y)T^ = 1 S(y) (t^ - 3M 4 2 D 4 7r) , (35) 



and so will remain finite on the codimension- two brane, but is plagued with a much worse pathology: 
7r has the wrong sign kinetic term on the brane in the UV. Indeed in that regime, the left-hand side 
of (|35|) is negligible, and ir satisfies tt = 1/(3M|D4)T, while in the IR the left-hand side dominates 
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and 7r exhibits a 6D behavior and couples to —T. Its kinetic term hence changes sign, signaling 
the appearence of a ghost. In analogy with the transverse, traceless case, its solution is 

vr=Al " -r (4) , (36) 

M 4 2 lD 4 + 5 (-n 4 ) 12 y 1 

and its propagator is thus not positive definite. Notice that for the tensor modes, the CU appears 
in (|34h with a negative sign while it comes in with a positive sign for the scalar mode (136 j) . 

The ghost manifests itself in the gravitational exchange amplitude A between two conserved 
source T^ u and T' on the brane: 

A ~ f d 4 X fe) TT + TT^) T"- 



A / d 4 x ' 1 



Mi J _n 4 + 5 (_n 4 ) 



5- " 3^ T ') + S in 4 + 1 ,(-n 4 )^ T ' } • (37) 

The tensor structure is recognized as the sum of a massive spin-2 contribution, with the famous 
1/3 coefficient, plus a conformally-coupled scalar. In the IR, corresponding to — Q4 <C g(— O4), 
the overall coefficient of the TT' term asymptotes to —1/4, as expected from the six-dimensional 
behavior of the force law. In the UV, however, the overall TT' coefficient tends to —1/2 which 
indicates a pathology. In that regime, the scalar amplitude is indeed negative, indicative of a ghost. 
The presence of this ghost is not specific to our cascading framework and appears in general 



codimension-two and 
on the 3-brane. In 
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righer scenarios 1CJ, [llj • To remove it we can modify the operators localised 
it was argued that adding a tension to the codimension-two brane removes 
the ghost providing the tension satisfies To > |M 4 mg. However, this tension is also necessarily 
bounded from above by To < 2ttMq = 2 r Kvnz > m§M\, due to the normal requirement that the deficit 
angle of a codimension-two object must be less that 2tt for the induced metric to remain flat. 

In what follows we shall deal with an alternative way to remove the ghost, which does not require 
the introduction of a tension on the brane, but achieves the same effect of adding additional kinetic 
terms for ir from brane localized kinetic operators. 

V. BRANE REGULARIZATIONS 

As mentioned in the previous section, the presence of the ghost ir is generic to any codimension- 
two and higher framework [lcl. [ill]. One can indeed recover the pure codimension-two scenario by 
considering the limit where the five-dimensional Planck mass vanishes. In that case mg 3> \/— D 4 
and the gravitational exchange amplitude is still given by (|37p . with 

,(-□2) = ^^^=) . (38) 

Here, A is a UV cutoff that now needs to be introduced since the geometry on the codimension-two 
brane is no longer regularized by the codimension-one brane. In other words, A represents the scale 
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at which new physics is involved in the description of the brane. Considering the codimension-one 
brane as a regulator, we would have A = 2Mg/M| . We first focus on the regularization of the pure 
codimension-ra scenario (no cascading), but the results are fairly general and will be extended to 
the cascading framework in the next section. 

These inevitable divergences in higher codimension force us to regularize the branes. There 
are many ways to go about this, ranging from the elaborate construction of smooth topological 
defect models of the brane, to crudely imposing cutoffs. In the gravitational case it is crucially 
important that any chosen regularization scheme be gauge invariant. In Appendix [Bl we propose 
two alternative regularization schemes that all lead to the same result, as well as giving a description 
of how we understand the absence of a ghost by Kaluza-Klein reduction within the brane. We review 
these in the rest of this section after first discussing how we can always construct gauge invariant 
regularizations by means of introducing Stiickelberg fields. 



A. Stiickelberg Approach 

A straightforward way to regulate branes is to replace them with 'smeared branes' or 'brane 
distributions'. Essentially we can think of a brane as being made up of N identical branes located 
at the same position. To smear them we spread the branes out so that we have some distribution. 
To understand this, let us first consider the more straightforward example of a codimension-one 
brane. Imagine the brane is located at <3?(x,y) = 0, where $ is some scalar function, then the 
action for the brane will be 

Sbrane = / <^X y/-giC M (V ' X*) ) ( 39 ) 

where x % denotes matter fields localized on the brane, and is the usual induced metric on 
the surface $ = 0. This action is invariant under the full diffeomorphism group since 5D dif- 
feomorphisms act as 4D diffeomorphisms on the surface of the brane. We now smear this into a 
distribution /(e) of branes localized at <J?(cc,y) = e, such that $™ def(e) = 1, and define g e ^ u as 
the associated induced metric. We thus have 

/oo r 
def(e) / d 4 x^£ M (^, xl • (40) 
-oo J 

In general there is a 4D tensor 7^(2,2/) such that 

9e,i*u = l l iu(x,y)\<s>= e , (41) 

so we may write 

/oo f f 

def(e) J d 4 x J dyd y $6($ - e)^C M (^ v , X l ) ■ (42) 

Performing the integral over e we obtain 

CLe = j d ' x j dydy^f^^Cu^x 1 ) ■ (43) 
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Since the original action was gauge invariant, then so is this regularized form. However in practice 
it is most simple to understand this in the gauge in which $ = y so that 

s ZLe = J [ d ' x J dyMV^giCvig^x*) • (44) 

Following the same steps for codimension-two branes, defined as the locus of 3>i = 0, $2 = 
we have 

CLe = J [d 4 x J dydzi^^f^u^V^CM^,^) , (45) 

where {$1, $2} = dy<&\d z $2 — d z $idy&2- Again this is most straightforward to understand in the 
gauge where <3?i = y, $2 = z so that 

5 brLe = J d ^ x J dydzf(y,z)^glC M (g f , u ,x t ), (46) 

or, in polar coordinates, 

CLe = J " d *x J rdrdOf^^lCM^xl • (47) 

Intuitively we can think of <I>i and $2 as Stiickelberg fields corresponding to the gauge transfor- 
mations which are spontaneously broken by the branes. In our calculations we mostly choose de 
Donder gauge in the bulk, which gives us sufficient gauge freedom to set the branes at a fixed 
position. This is equivalent to working in the gauge in which the Stiickelberg fields take the values 
<&i = y and $2 = z. 

Having established the gauge invariance of the smoothing procedure, in the rest of this section 
we present two alternatives regularization schemes that lead to the same low energy physics. 



B. Spherical Regularization 

In the first scheme, we replace the codimension-n brane with a codimension-one brane, wrapped 
around (n — 1) compact directions (see [13I ] for a very similar construction). Expressing the 
Minkowski background in terms of the generalized spherical coordinates along n directions, 

ds 2 D = dr 2 + r 2 d^_ i + d ^ Xi ^ ( 48 ) 

with D = 4 + n, and r = corresponds to a codimension-n point. If, instead, we choose to localize 
the brane at r = A, we are then dealing with a codimension-one brane. And since this brane is 
now (3 + n)-dimensional, the corresponding intrinsic Einstein-Hilbert is then Rs+ n rather than R4. 
As we can see in the Appendix, this has dramatic consequences for the scalar degrees of freedom. 
In particular, the induced metric perturbations sourced by an energy-momentum T^ v on the brane 
is now given by 

V = (D 4 ) ( V - ^— T7] \ + . . . (49) 
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where the ellipses indicate total derivative, and the propagator Ga^a) tends to ~ — I/M4CI4 in 
the UV and to (V - O4/ A) 2 ~ n /M^ in the IR. It should be clear from this result that the resulting 
gravitational amplitude is ghost-free. Notice, furthermore, that we obtain the same coefficient for 
the TT 1 amplitude at any energy scale. Strong coupling effects will however modify this coefficient 
at high energy, as we discuss in Sec. I Villi 

C. Medium Model 

In the second scheme, the codimension-n brane is promoted to a full D = 4 + n-dimensional 
object whose thickness A now plays the role of the regularization scale. To realize this, we divide 
the entire space-time into the two half-spaces, r < A and r > A, where the in-space models the 
blurred codimension-n brane. 

Since the brane is now D-dimensional, the natural curvature invariant to include on its world- 
volume is Rd instead of -R4. In all generality this can have a different Planck scale from that of 
the ambient space-time. In other words, the brane is thought of as a medium with gravitational 
permeability e = M 2 /M® -2 [15]. The gravitational part of the action is therefore given by 

S = ^— d D X R D + g / d D Xy/=gER D . (50) 

^ Jr>A ^ Jr<A 

In this case the brane localized R4 term arises naturally from the D-dimensional Ricci curvature 
in the thin-brane limit A — ► 0. 

Although the technical details in this case are slightly different, as shown in the Appendix we 
still recover the same expression for the brane induced metric perturbations as (|49p . giving rise to 
the same ghost-free amplitude. 

D. Effective Theory Approach 

The spherical regularization and medium approach discussed above offer two independent 
smoothing procedures that lead to identical ghost-free amplitudes. Here we check this idea by 
deriving a low-energy effective action in the thin-brane limit, obtained by doing a KK reduction 
over the world- volume of the regularized brane. 

Indeed we find that either choice leads to the same effective action as the thickness is taken to 
zero. Let us therefore consider the medium action (|50[) for concreteness. Within r < A, we take a 
metric ansatz (we can ignore the KK vectors since they are not sourced) 

ds 2 r<A = g^dafda? + e 2<t> (dr 2 + r 2 df^„ 5 ) , (51) 

where g^ u and are functions of x^ only. In this approximation, we have Rd = R4 — 2{D — 
4p 4 </> - (D - 4)(D - 3)(d^4>) 2 . Substituting into and integrating over the radial and angular 
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directions, we obtain, after integration by parts, the following effective action 

M D ~ 2 (' /If 2 f 

S = — ^— / d D xV=gER D + ^ / d 4 x^e n * (i? 4 + (D-A)(D- 5)(<90) 2 ) , (52) 

1 J Af out 1 J 

where M| ~ eMj^~ 2 A n is the induced 4D Planck scale. As mentioned earlier, exactly the same 
effective action follows from the spherical regularization. 

The effective brane action therefore reduces to an DGP term, although with additional 
scalar-tensor couplings to the extra-dimensional modulus 4>. Note that in the standard DGP case, 
corresponding to D = 5, the kinetic term for (f> vanishes, while the extra factor of in the measure 
can be set to unity by working in Gaussian normal coordinates. In other words, the modification 
is pure gauge in this case, leaving the standard DGP framework unaltered by the regularization. 

However, the scalar-tensor structure is crucial for higher-codimension branes — the extra con- 
formal coupling affects the TT' part of the exchange amplitude, thereby curing the ghost. As 
shown in detail in the Appendix, by following similar steps to those of Sec. HVI we arrive at the ten- 
sor structure displayed in (|49p . with the TT' coefficient appropriate for Z)-dimensional embedding 
space-time. In what follows we apply these regularization techniques to the cascading framework. 

VI. GHOST-FREE CASCADING GRAVITY 

To proceed we apply the general approach of Sec. |V] — first smoothing out the brane to finite 
thickness and identifying the correct low energy theory. The resulting action has a four-dimensional 
DGP term, albeit with additional scalar-tensor couplings involving the extra-dimensional metric 
components. For concreteness, we henceforth focus on the codimension-two set-up; the generaliza- 
tion to higher codimension is straightforward. 

Remarkably, these new scalar-tensor couplings eradicate the ghost, leaving an exchange ampli- 
tude that is manifestly free of instabilities. The corresponding force law extrapolates between a 
six-dimensional (~ 1/r 4 ) behavior in the IR to a four-dimensional (~ 1/r 2 ) behavior in the UV. 
Meanwhile the tensor structure is that of six-dimensional gravity on all scales. 

Thus we begin by smoothing out the codimension-two brane to have finite thickness. As argued 
in Sec. El it no longer make sense to add an term since the worldvolume is now 5+1 dimensional. 
Instead, the appropriate brane-induced gravity terms are either an R§ term on the (codimension- 
one) boundary of the regularized brane, as in Sec. IV Bl or an Rq term throughout its worldvolume, 
as in Sec. IV Gl 

As argued in Sec. IV T)\ either choice leads to the same effective action (|52p as the thickness is 
taken to zero, which for D = 6 reduces to 

Seff =^Y j d^x^ie 2 * [R 4 + 2{d4>?) + S mattCT [g] ■ (53) 
Note that the matter action is insensitive to the regularization, as argued in Sec. IV Al Thus the 
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effect of regularization amounts to modifying the brane gravitational action to a scalar -tensor 
theory. 

Conservation of stress energy requires that the variation of this action with respect to <p vanishes, 
which, at the linearized level, gives the condition 

i? 4 = 2D 4 . (54) 

Once again we can expand the 4D metric into scalar and TT components as in (j32l) . The TT part, 
hjjj , is by definition oblivious to the trace part of the brane source, and is hence unaffected by the 
modified form of the brane effective action. The scalar-tensor structure does, however, affect the 



trace mode it. Indeed, the trace of the effective source from ()53[) is 

T eff = 6(y) [m 2 4 {R a - 6D 4 <£) + T^} = <%) {-2M 4 2 i? 4 + T< 4 »} . (55) 
Furthermore, in terms of the decomposition (|32p we have i? 4 = — 3D 4 7r. Hence, (|35p becomes 



2 



1 „, M 



2 



□ 5 - mev^SsJ vr = — 5{y)T + ^-5{y)U^ . (56) 



Note that the kinetic operator now coincides with that in the TT equation. In particular, the 
brane-induced kinetic term for ir is healthy, which traces back to the additional [ZI 4 </> term in (|55f) 
flipping the sign of i? 4 . The solution for tt, 

2 1 1 



7T 



• -r< 4 ) , (57) 



Mi -D 4 + gi-Di) 12' 
is therefore exactly of the same form as its TT cousin. 

Thus the one-particle exchange amplitude between two conserved sources and T' on the 
codimension-two brane takes the simple form 

A ~ M I ** T "' -□,/,(-□.) " M ■ (58) 

displaying the tensor structure of six-dimensional gravity. This is manifestly ghost-free since the 
coefficient of the TT' part satisfies —1/4 = —1/3 + 1/12, corresponding to a massive spin-2 plus 
a healthy conformally-coupled scalar. Notice that the above tensor structure is such that a cos- 
mological constant does not gravitate at the linearized level. This is consistent with the fact that 
pure tension on a codimension-two brane leaves the induced geometry flat and generates instead 
a deficit angle in the embedding space. Note, however, that this codimension-two mechanism for 
taming the backreaction of vacuum energy is different than degravitation. 

This tensor structure is of course different from the usual 1/2 predicted by general relativity. 
However, as in standard DGP, we expect that non-linearities decouple the extra scalar modes in 
the vicinity of astrophysical sources, so that the theory approximates Einstein gravity locally. In 
particular, in the limit rag < m§ in which gravity has a 5D regime at intermediate distances, we 
expect a double Vainshtein effect: as we approach a source from infinity, first n then hJ^J become 
strongly coupled, while the TT' coefficient goes successively from 1/4 to 1/3 to 1/2. We will come 
back to this issue in Sec. IVIIII 
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VII. CURING THE GHOST IN THE DECOUPLING LIMIT 

It is enlightening to study the origin of the ghost in the decoupling limit of the theory, which 
corresponds to sending M^,Mq — ► oo while keeping the strong coupling scale (Mg 6 /M|) 1//7 fixed. 
The bulk physics can be integrated out in this limit, resulting in an effective 5D description in 
terms of a metric h a p and a scalar mode tt encoding brane bending in the extra dimension. 

At quadratic order the 5D Lagrangian is given by 

£s = {~ ha(3 ( £h U ~ 3vr (-n 5 /i 7 7 + d a dph a ^ } + 5{y)U , (59) 

where h a p is the induced metric on the co dimension-one brane, and £4 describes the codimension- 
two brane. Since tt is a scalar, £4 can in principle depend both on the induced metric h^ u , as well 
as on this scalar mode. Thus let us assume the general form 

ft/r 2 1 

£4 = -f {-h^(£h)^ - (3tt {-Dih^ + d^h^) + ZynU^} + , (60) 

where (3 and 7 are constants. In particular, the naive DGP extension studied in Sec. IIVI corresponds 
to p = 7 = 0. 

Following Appendix A, we find that tt appears on-shell in the decomposition of the 4D metric 

as 

hfja, = hjj + nr]^ + gauge terms . (61) 



This agrees with the earlier decomposition — e.g. (|32p — , hence the choice of notation. Moreover, 
it immediately follows that tt is related on-shell to the induced scalar curvature by R4 = — 3D47T. 
Using this fact, the equation of motion for tt is thus 

]\/[3 1 ]\^2 

- -f^TT = -6(y)T® + -fS(y) (0 + 7-1) ^ , (62) 
from which we infer the following structure for the one-particle exchange amplitude 
rpuv 1 / , 1 , 1 -D 4 + m 5 v / -D4 

-□ 4 + m 5V / -n4 \ 3 12 _q 4 ip+i 11 + m5 ^/rn^ j 

Therefore, the coefficient of the trace part interpolates between —1/4 at the largest distances, 
corresponding to the tensor structure of 6D gravity, and 

-K'-^n) <«j 

at the shortest distances. 

For instance, in the case f3 = 7 = of Sec. IIVI the above equation for tt agrees with (j35|) in 
the me — ¥ limit — the induced wave operator has the wrong sign, indicative of the ghost. More 
generally, though, the tt mode has a healthy propagator whenever 

/3 + 7>l, (65) 
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which therefore requires non-zero (5 and/or 7. 

Let us therefore revisit our regularization schemes under this light, in particular in the effective 
theory language of Sec. IVI1 Since R4 = — 3D47T on-shell, comparison with (|54p implies that 2(f) = 
—3tt. Substituting this in the expansion of (|53p to quadratic order, we find that our regularization 
prescriptions all correspond to f3 = 6 and 7 = — 3. In this case ([621) reduces to 

- ~fa 5 n = -5{y)T& + -f5{y)U^ , (66) 



again in agreement with the m§ — > limit of (|56j) . 

Thus we see that our different ways of smoothing the codimension-two brane are just one 
example of how the ghost can be cured. Of course not every choice of (3 and 7 can necessarily 
be embedded within a covariant 6D theory. Another example where this is possible, however, is 
through the addition of a tension A on the codimension-two brane 12| . Thanks to tt self-interactions 
in the non-linear extension of (|59p . the non-zero tension generates an effective kinetic term for tt 
localized on the brane, corresponding to 

1 = l^M~ (67) 

And since f3 = in this case, ([65]) yields a lower bound on the required tension. 

VIII. STRONG COUPLING PHENOMENON 



Since massive gravitons in Minkowski spacetime exhibit the famous vDVZ discontinuity 171 ]. 
due to the persistence of an additional scalar degree of freedom in the m — > limit, any theory in 
which the graviton is a resonance must exhibit strong coupling of any additional scalar modes at 
some scale, known as the Vainshtein scale, so as not to be ruled out by standard tests of General 
Relativity. In the 6D cascading model this means that the factor of —1/4 must effectively be 
converted to —1/2 in the strong coupling region. Here there are two additional scalar degrees of 
freedom, which intuitively arise from the h zz and h yy components of the metric. The scalar arising 
from the 6D — > 5D transition is tt, and we expect it to become strongly coupled according to the 
natural generalization of the standard DGP story generalized to 6D. The scale at which strong 
coupling occurs, i.e. the Vainshtein scale, can be determined from the equation of motion for tt in 
the decoupling limit, including the nonlinear terms for tt 

~ DsTT + ^| {{d a d^f - (Dvr) 2 ) = 8(y)-±g (t + M^) , (68) 

by comparing the magnitude of the terms nonlinear in n, with the linear ones. To determine this 
we may first solve the linear equation, whose formal solution for y > is 

tt(x^, y) = e- v/rni %o(^) • (69) 
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Here ttq is determined from the boundary condition imposed by the delta function source, d y ir = 
-5BF( r + M ? fl(4| ).t°be 

Since i?( 4 ) contains terms of order D^n, there are two regimes of interest. For distances r > «5 
we have ttq ~ ^g^-j T, whereas for r <C m?, we expect ttq ~ M \^ T. 

From the equation of motion, we see that strong coupling kicks in when ddir ~ {ddn) 2 /m|. In 
terms of the source this is when 

r > m 5 dT ~ m§Mf , (71) 
r < m 5 T ~ m^M| , (72) 

and so for a mass M with T ~ Mj r 3 , the Vainshtein scale is given by 



r>>m5 ^ = V^M|' (?3) 



r«m 5 Rv = ^^- (74) 

Once the n mode has become strongly coupled, the effective 5D theory reduces to pure gravity 
coupled to a codimension-one brane with induced kinetic terms. In other words, it looks like the 
standard DGP model, and so we anticipate that the second additional scalar will become strongly 
coupled at the usual scale 



*" = *VSP3- (75) 

For this argument to be consistent we require Ry > Rv- For distances r <C ms, this is automatically 
true provided mg < 7715 — a necessary condition for a cascade from 6D to 4D via 5D. 

Although the details of the strong coupling transitions should be checked by a careful calculation 
of the effective two field 7r-Lagrangian, we can see at least qualitatively that there can be a well 
defined cascade of strong coupling scales provided tuq < m^. We should stress that there is no 
reason to doubt that the same mechanism will also occur for m§> m^. However, it will be necessary 
to do more explicit calculations in this case as there may be some subtle interaction between the 
two scalar modes which is not accounted for in the above argument. 
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APPENDIX A: CASCADING GRAVITY A NAIVE APPROACH 



Starting from the 6D action (|31l) . and considering perturbations around a flat metric 

ds 2 = (tjab + h AB )dx A dx B , (Al) 
in de Donder gauge <9a^b = \dsh^ the 6D Einstein equations take the simple form 

MlGf B = -M De (h AB _ he ^ = _ 6{z) ( M f G§ _ T (5)) s * a5 P b . (A2) 

Here the five-dimensional stress-tensor, T^) , includes the four-dimensional Einstein term, G$ , as 
well as matter stress-energy localized on the (3 + l)-brane, T^): 



r (5) 
■a/3 



6(y) (-M 4 2 GW + Tj0) 6^ . (A3) 



1. Step 1: 6D -> 5D 

Since there are no energy-momentum sources along the z direction, the (a, z) and (z, z) Einstein 
equations allow us to set h az = and h zz = h^, respectively. The induced five-dimensional gauge 
choice is thus d a h°p — dph a a = 0, in terms of which the five-dimensional Einstein tensor takes the 
form 

G S = ~\ ( D 5 Kp - dadphfy . (A4) 
The remaining five-dimensional Einstein equations are then 



Mi 



□ 6 + ® Hkj h aP = 5(z) - V J\ - ^-Mid a dph% , (A5) 



2 

with ttiq defined in (|19[) . In particular, the trace part satisfies 

Ma 1 

--fUM, = — A S(z)T^. (A6) 

Note that, since ~ <5(y)> this mode only couples to the co dimension- two source. That is, 
it does not "feel" the codimension-one brane, and its solution therefore suffers from logarithmic 
divergences characteristic of codimension-two branes. However this is purely a gauge artifact as 
can be seen by going to the gauge defined by h z M = 0, where the metric is given by 

~ K e = h^f TT + J£> %^T( 5 ) . (A7) 

«P a/3 -r ^ V ) 

In the above we have found it convenient to project out the 5D transverse and traceless variable, 

,(5d)TT , dad/j 



24 



which satisfies the decoupled and fully regularized equation 



(A9) 



This can be solved exactly following the 6D — > 5D reduction in the scalar toy model — see QUI Aft : 



Mi 



□5 - m 6 y/-\J 5 



h 



(5d)TT 

a/3 



T< 4 > - M 4 2 G (4) 



■4 I ^ 



T< 4 ) + M|i?W 



(A10) 
S(y) , (All) 



where h^? TT now denotes the metric fluctuation evaluated on the codimension-one brane, and 



where we have substituted for T^) using (|A3p , Note that the derivative operator on the left-hand 
side is recognized as the inverse of the scalar Green's function (see Sec. IIII A|) . 



2. Step 2: 5D -> AD 



Writings = -^hh> yy 



□ 5 ,(6d)TT _ □. 



3m' 



r-rj^h^u T , we can express the 4D transverse and traceless part 



as 



^(4d)TT = h (5d)TT 



3 §ir 4 ) k^ dv ^ 



-15 / LJ4 

giving the following expression for the four-dimensional Einstein tensor 



□4 



7T 



Meanwhile, the (y, y) component of (lAlOj) yields an equation of motion for ir: 



M 3 



1 



□5 - mev^J vr = — <%) [T^ - 3M 4 2 n 4 7r] . 



(A12) 



(A13) 



(A14) 



Putting everything together, the {j-iv) component of (|A10|) simplifies to 



M 3 r 



Ml 



□5 - m 6 v^J h^ TT = 8{y) ^ + ^ D 4 , (A15) 

where = — ^T^rj^ + ^j^T^ 4 ) is the transverse and traceless part of the matter stress 
tensor. 

Equation (|A15|) is the exact analogue of the scalar field equation derived from (fl~3j) . Thus, as 
advocated, the transverse and traceless story parallels the scalar toy model presented earlier. The 
solution is therefore given by 

2 1 



h (4d)TT 



Mi-D A + g{-a 4 ) 



J fJ,U ) 



(A16) 



where the function g was introduced in ([2"3"j) , and h^v now denotes the induced metric fluctuation 
on the codimension-two brane. 
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Equation (|A14p for the scalar mode n, on the other hand, is slightly different as that for the 
tensor modes and gives rise to the solution 

1 1 



71 ~ 3M|-D 4 -2 5 (-D 4 ) T(4) - (A17) 
A key difference between (|A16p and ()A17p is the presence of the negative sign for g{— O^) in the 
expression of w, giving rise to the presence of a ghost as soon as D 4 S> g{— CI4). 

APPENDIX B: CODIMENSION-n BRANE REGULARIZATION 
1. Spherical Regularization 

Let us work in D = (4 + n) dimensions where n is the codimension. Consider the action 
representing the codimension-n version of the DGP model: 

S = -EL— J d D x^g-ER D + ^ L J d A x^g~4R4 + J d^x^id m . (Bl) 

In this section we replace the codimension-n brane with a codimension-one brane which is wrapped 
around (n — 1) compact directions (see [3]). 

In particular, let us assume that the bulk metric is given by Minkowski spacetime. It is useful 
to write this in the generalization of spherical or polar coordinates in the n directions as 

ds 2 D = dr 2 + r 2 d^_! + d 2 x A , (B2) 

where by d£l 2 l _ l we understand the line element on a unit (n — l)-sphere. The point r = is 
a codimension-n point where we assume the brane to be located. To regulate we replace the 
codimension-n brane with a codimension-one brane located at r = A so that the induced metric 
on the codimension-one brane is 

induced = A 2 d^_! + d 2 X 4 . (B3) 

In this way we can take our matter fields to be spread out over the codimension-one brane since 
in the limit A -> we are essentially performing a KK compactification and the dynamics on 
the brane will be dominated by the KK zero modes which are effectively living in D dimensions. 
The KK reduction approach will be pursued in Sec. IB 31 In particular this is true for energy on 
the brane E <C A -1 . Strictly speaking we should address the stability of putting the brane here, 
however it is likely easy to cook up a stabilization mechanism using e.g. form fields. Nevertheless 
we can still address whether this set-up has ghosts independent of this, since even an unstable 
solution should be unitary. 

Defining 7 -1 = A n_1 V^_i where V n -\ is the volume of the (n — l)-sphere, the regulated action 

is 

S = J d D xV^R D + J ^f J dP^Xif^lRD-! + 7 J d D - l x^g~^r x C M • (B4) 
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In practice we can neglect the KK excitations of the matter fields. For n > 3 we include a tension 
in the spherical directions to support the background solution for otherwise the extrinsic curvature 
will pick up a jump due to the nonzero curvature in the spherical directions. Since this is pure 
tension, however, it will not contribute in the perturbed equations, and its only role is to support 
the background solution. Furthermore, the tension is needed for regularization purposes only and 
disappears in the thin-brane limit. 

The perturbed equations of motion are thus given by 

mE~ 2 G (D) A b + S(r - A) 7 M 4 2 G (D ~ 1) \ 5 A J b B = 7 5(r - A)T% 5 A J b B , (B5) 

where A, B, . . . span as usual over all D dimensions, whereas a spans over all but the transverse co- 
ordinate r. It should be evident that this is invariant under the linearized version of D-dimensional 
diffeomorphisms. 

The reduction from D dimensions to (D— 1) is now familiar. We again choose the D-dimensional 
de Donder gauge 

V A {h A B - \5 A B h c c ) = . And since there are no sources in the (r, A) directions 
we infer as previously, h r a = and h r r — \hP c = , which implies h r r . = h a a . The de Donder gauge 
condition is then V a (h a b — 5 a b h c c ) = 0, which implies ^ = — \ (Do-ih a b — d a dbh c c ) . Since this 
is automatically traceless we infer that 

- \m%-\2 - D)U D h c c = 7 5(r - A)T , (B6) 

and so by defining the traceless part of h a b as 

H ab = h ab - -^—V a V b h c c , (B7) 

we have, after substituting in the trace, 

- ±M°- 2 n D H ab - i 7 «5(r - A)Mia D ^H ab = jS(r - A) (r ab - ^L_n ab T^ +..., (B8) 

where the ellipses indicate a total derivative, which are irrelevant when coupling to a conserved 
source. In the limit A — > we have Dd-i —> according to the conventional KK prescription. 
Defining the Green's function Q as the solution of 

Mg- 2 n D g(x; x') + 7 5(r - A)M 4 2 D^Or; x') = 5 {D) (x; x') , (B9) 

we infer that the gravitational exchange amplitude is given by 

A ~ j d 4 xh^T'^ ~ J d 4 xH^ v T'^ 

~ Jd 4 xg A (x^! [ T^--^_r l ^T^, (BIO) 
where £?a(^ m ) = G(x^,r = A; 0, r' = A). This is manifestly ghost free since we can write it as 

A ~ J d 4 x g A ( x ») | (t^ - ± Vijh ,t) t'^ + err' | , (bh) 

where C = \- ^ > for D > 5. 
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2. Medium Model for a Higher-Codimension Brane 

In this regularization, we consider the two half-spaces 

M- m = {X | r < A} , M out = {X | r > A} , (B12) 

with r denoting the Euclidean distance in the transverse direction and A the thickness of the 
in-space in the transverse directions. The in-space models the blurred codimension-n brane. We 
might think of it as a medium with gravitational permeability e = M 2 /M^~ 2 . The action for this 
set-up is thus 

S=-^— d D X^g^R D + D / d D Xyf^ER D , (B13) 

1 J Mont 1 JMin 

which corresponds to integrating over a sharp profile P(r) = eO(A — r) + 9{r — A) that describes 
the support of y 7 — <?d Rd in M\ n and M ou t- 

Varying the above action gives rise to Einstein's equations in both half-spaces, Mj n and M out , 
along with the following boundary conditions 

e n A J A \ 9M = n A J A \au • (B14) 

These are the analogue of the Israel jump conditions. In the weak-field approximation, the current 
density through the interior surface becomes J A = d A hg—dBh AB . Imposing the 4-gauge conditions 
d a h a = 0, the Einstein tensor in the (4 + n)-split is given by 

{8 D h)^ v = - - ((£ 4 /i) Mi , + (dfj,d u - t^CU) h n 

+n n [h llI/ - V ^(h 4 + h n )}+ V ^d a d b h ab ^ , (B15) 

{£ D h)ab = ((°4 + Dn) Kb + [d a db - Vab (p4 + Ql4 + K) 

-d c d ia h b)c + Vab d c d e h ce - d»d [a h b)fl + Vab d^d u v) > ( B16 ) 

where the traces are defined as follows /14 = r}^ v h^ v and /i n = r] ab h ab . The mixed components of 
the Einstein tensor are not required. 

We are interested in the polarization tensor of the induced weak field h^ u . In the case of induced 
gravity, the integration theory of the weak field equations of motion requires us to first gauge shift 
to a new field ^> ^ u that obeys a generalized Fierz-Pauli equation after appropriate gauge fixing. 
The form of the gauge-shift is clear from (|B15j) : 

hfj, v = + Ch n r]^ u , (B17) 

with C denoting an adjustable constant. Imposing further the D — 4 gauge conditions d a h a b = 
-d b h n /2, we find for C = -1/2 

(£4*)^ - m\U n ) (V^ - ¥4^) = ^=2^ , (B18) 
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where we have allowed for an (induced) source term. This is indeed Fierz-Pauli theory for the 
shifted gauge field with a generalized mass term m 2 = — □„ in accordance with the usual 
Kaluza-Klein reduction. 

Since only transverse degrees of freedom propagate in Fierz-Pauli theory, the trace of the gauge 
field can be eliminated for a conserved source, ^ 4 = 2T 4 /3M^~ 2 D n , allowing us to solve ([B18]) : 

= -j^ n^mo ( T - - \^ T ) ■ (B19) 

The boundary conditions (|F314|) impose m 2 (D n ) = /(CU). However, \I> ^ depends on h n . In order 
to eliminate h n we write down the equations of motion for h a b- 

(□4 - m\U n )) (hat + IhnVa^j ~ d»d {a h b )„ = ^|zf (vab ~ =J^j ) T ' < B20 ) 

Tracing (1B2U1) and using again d a h a ^ = gives h n entirely in terms of the induced source term 
and the appropriate Green-function: 

4 D-5 T 4 
n ~ 3M^ D^2 U 4 -m 2 {U n ) ' (B21) 

Note that for D = 5 the source term vanishes. For D > 5 the solution (|B21|) guarantees that hu V 
has the correct polarization tensor. Indeed, 

hilV = "M£- 2 D4-m 2 (n n ) " D^2^ uT ) • (B22) 

3. Dimensional Reduction in Zero thickness limit 

As argued in Sec. IVJ both regularization schemes described above lead to the same 4D effective 
brane action in the thin-brane limit. To review, we perform a KK reduction within the world- 
volume of the thick brane by taking the ansatz ds 2 D = g^dx^dx" + e 2 ^ [dr 2 + r 2 dQ 2 j_ 5 ) . The 
resulting effective action is given by 

M D ~ 2 r m 2 r 

S = — |— J d D x^9BR D + d A xy/=g^e n * [R 4 + (D - 4)(D - 5)(S0) 2 ) . (B23) 

A key difference with the usual DGP Einstein-Hilbert term on the defect are the extra terms 
involving <f>. These will prove instrumental in getting the desired tensor structure for the exchange 
amplitude. 

As previously, in de Donder gauge the Einstein equations take the form 



m£- 2 /. , 1 



, n D {h A B -~6 A B h c c )=T A B ^. (B24) 

The 0-dependent terms in the brane action contribute to the (r, r) component of the stress tensor. 
Conservation of energy-momentum therefore requires this contribution to vanish, which, at the 
linearized level, yields the condition 

R 4 = 2(D - 5)D 4 </>. (B25) 



Meanwhile, the stress tensor components along the brane are given by 

T A* (6) = tf (n)( r ) |_ M | G ^(4) _ p _ 4 ) M 2 _ 3^ ^ + T£ j . 

Now the choice of gauge in the bulk is consistent with the gauge relations d u h 
h^jj, = 2(6 — D)(p. In particular, the linearized 4D Einstein tensor takes the form 

G$ = ~\ (D 4 V + " 5 P4^ - 2(Z> - 4)«9 M cU) , 

whose trace is indeed consistent with (jB25]). 

Substituting into the (fx, v) equations of motion, we find 

- 1 (m£- 2 D 6 + «5W(r)M 4 2 D 4 ) {K - 28t<f>) = 8^{r)Tt . 
The scalar term can be eliminated by taking the trace, 

- (m£- 2 D 6 + 5(")(r)M|D 4 ) 4> = -^Z^5 (n) (r)T, 
and substituting the result back into (|B28P to obtain 

- \ (m^- 2 D 6 + ^)(r)M 4 2 D 4 ) K = 6^(r) - q^KT) • 
This agrees with our earlier results. 
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